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\ Abstract. We constrain the number density and evolution of Compton-thick Active Galactic Nuclei 

\^ . (AGN), and their contribution to the extragalactic X-ray background. In the local Universe we 

use the wide area surveys from the Swift and INTEGRAL satellites, while for high redshifts we 
explore candidate selections based on mid-IR parameters. We present the properties of a sample 
of 211 heavily-obscured AGN candidates in the Extended Chandra Deep Field-South (ECDF-S) 
selecting objects with /24|j„,//i;> 1000 and7?-/r>4.5.TheX-ray to mid-IR ratios for these sources are 
. significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for 

the general AGN population, suggesting column densities of Nh>5 x 10^^ cm^^. The space density 
of CT AGN at z^2 derived from these observations is ^lO^^Mpc^^, finding a strong evolution in 
the number of Lx>lO'^'^ erg/s sources from z=1.5 to 2.5. 

2 ■ Keywords: galaxies: active, Seyfert. X-rays: diffuse background 

^ ■ PACS: 98.54.Cm, 98.70. Vc, 98.62.Js 



O 



(N 



: COMPTON THICK AGN WITH INTEGRAL AND SWIFT 

> , 

The most obscured AGN known are those in which the neutral hydrogen column 
[-^ '. density (Nh) in the line of sight is higher than the inverse Thomson cross section, 

CN ; Nh^I-SxIO^"^ cm"2. These are the so-called Compton-thick (CT) AGN. Contrary to 

CN ■ the situation for less obscured sources, not much is known about the number density 

^ . of CT AGN. Thanks to the deep Chandra and XMM-Newton surveys it is now clear 

I that the fraction of moderately obscured, Compton-thin, AGN is on average ~3/4 of all 

AGN, and is higher at lower luminosities [1, 2] and higher redshifts dH], but there are 
no comparable constraints on the number of CT AGN. One of the best ways to find CT 
I AGN is by observing at high energies, namely E>10 keV. One clear advantage of these 

o3 I high-energy observations is that photoelectric absorption has minimal effects, so even 

CT AGN can be easily detected. 

Using the IBIS coded-mask telescope, INTEGRAL surveyed ~80% of the sky down 
to a flux of 5 mCrab in the 17-60 keV band. The catalog of Krivonos et al. [7] reports the 
properties of 130 sources detected in these all-sky observations and classified as AGN. 
Five of the 130 AGN are CT AGN. Similarly, Tueller et al. [5] presented a catalog of 103 
AGN detected in an all-sky survey with the Swi/i/BAr telescope. Excluding blazars, BL 
Lac and low galactic latitude observations, we obtained a sample of 89 sources, where 
only one source remains unidentified. In the Tueller et al. [5] catalog there are five AGN 
with estimated Nh greater than 10^^ cm^^. Figure [H shows the cumulative number counts 
of AGN, with CT sources shown separately, as a function of hard X-ray flux. In order 



to avoid the necessity of specifying a standard spectrum to convert fluxes to different 
energy bands, we show the INTEGRAL and Swift sources separately, but note that a 
good agreement (within ~40%) in the normalization between the two distributions exist 
if a standard band conversion is assumed. At these high fluxes the slope of the log N- 
log S is Euclidean, implying an uniform spatial distribution, as expected given the low 
redshifts of these sources. We also compare with the distribution predicted by the AGN 
population synthesis model with which Treister and Urry [0] fit the XRB, and find in 



FIGURE 1. LogN-logS distribution for 
AGN detected at high energies. The red line 
shows the AGN in the well-defined Swift/BAT 
samples in the 14-195 keVband [5], while the 
bottom panel shows the INTEGRAL sources 
101 in the 17-60 keV band. Solid squares 
show the 14 sources detected in the deep 3 
Msec INTEGRAL observations of the XMM- 
LSS field (S. Virani in prep.). Solid cir- 
cles mark the CT AGN detected with Swift 
{top panel) and INTEGRAL (bottom panel). 
The black solid lines shows the expected 
AGN logN-logS from the population synthe- 
sis model of Treister and Urry which at 
these fluxes corresponds to a Euclidean distri- 
bution. The dashed lines mark the Euclidean 
slope normalized to the number of Swift and 
INTEGRAL CT AGN. 



Since now the number density of CT AGN can be constrained independently, we can 
attempt to match the observed spectrum and intensity of the X-ray background (XRB). 
In Figure [21 we show our new fit, which matches the both the INTEGRAL and Swift/BAT 
observations at E>10 keV and the Chandra measurements at lower energies (which are 
~30% higher than the HEAO-l A2 observations). These new data confirmed that the 
original HEAO-l normalization should be increased by ~30% and ~10% at low and 
high energies respectively. In contrast, the AGN population synthesis model of Gilli 
et al. ^ assumed the original HEAO-l intensity at all energies, which translates into a 
relatively lower contribution from unobscured sources. In order to produce the necessary 
hard spectrum, Gilli et al. [6] had to assume a relatively high number of obscured sources 
at high luminosities, i.e., an unusual, inverted dependence of the obscured fraction of 
AGN as a function of luminosity [7J. 

IR SELECTION OF HIGH-z CT AGN 

Because most of the radiation absorbed at X-ray and UV wavelengths is then re-emitted 
in the mid- and far-IR, recent studies at these energies, mostly taking advantage of 



general good agreement in slope and normalization. 
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FIGURE 2. Observed spectrum of the ex- 
tragalactic X-ray background from HEAO-l 



H, Chandraj9],XMM 1 10], INTEGRAL HI 
and Swift U2R data. The dashed gray line 
shows the XRB spectrum from the AGN pop- 
ulation synthesis model of Treister and Urry 
101, which assumed a 40% higher value for 
the HEAO-l XRB normalization. The thick 
black solid line shows our new population 
synthesis model for the XRB spectrum; the 
only change is the number of CT AGN, which 
is reduced by a factor of 4 relative to the 
number in Treister and Urry |2]. Red, blue 
and thin black solid lines show the contribu- 
tion to this model from unobscured, obscured 
Compton thin and CT AGN respectively. 
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Sp/fzer observations, have been very successful in finding heavily-obscured AGN candi- 
dates missed by X-ray selection up to high redshifts. In particular, Fiore et al. iITbIi pre- 
sented a selection method based on the 24 /im to R band flux ratio and R-K color; specif- 
ically, /24//^>1000 and R-K>A.5 (Vega). According to simulations based on stacking 
of the X-ray signal, they estimate the fraction of heavily-obscured AGN in this sample 
to be greater than 80%. We applied this selection criteria to the MIPS -selected sources 
in the ECDF-S to select CT AGN candidates. Of the 7201 24/im sources detected to a 
flux limit of ~35 /iJy, 211 (~3%) satisfy the A4//i?>1000 and R-K>A.5 cuts. Of the 
651 X-ray detected sources in the Virani et al. [| 14| ] catalog, 18 are found in this region, 
~2.8%, a similar fraction as in the general population. 

As can be seen in Fig. [3l the X-ray to mid-IR ratio for the IR-red excess sources is 
about two orders of magnitude smaller than the average value for the X-ray-detected 
sample and for most sources falls in the range expected for obscuration of ~5x 10^^ to 
cm . The observed ratio for the IR-red excess sample is significantly larger than 
for the sources outside our selection region, even at similar rest-frame 12 /im luminosi- 
ties. The observedX-ray to mid-IR ratio for IR-red excess sources atLi2^OT~10^5 erg s~^ 
is roughly 3-4 times larger than the values expected from the X-ray versus mid-IR lumi- 
nosity for star-forming galaxies, thus confirming the AGN nature for the vast majority of 
our sample. The space density of CT AGN as a function of redshift inferred from these 
observations is shown in Fig.[3l 
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FIGURE 3. Left panel: Ratio of hard X-ray to mid-IR versus rest-frame 12 /xm luminosity. Red and 
blue circles show the obscured and unobscured X-ray detected sources respectively. Black triangles with 
error bars show the average ratio for these sources. Circles enclosed by squares identify the X-ray detected 
IR-red excess sources. The dashed line at a ratio of ^0.6 shows the average value of LxlL\2 for all 
the X-ray sources. The dotted line shows the relation between intrinsic X-ray and 12 jj.m luminosity 
for local Seyfert galaxies. Green circles show the location of the most luminous rest-frame 12 jim X- 
ray undetected sources outside the selection region. X-ray undetected IR-red excess sources are shown by 
olive circles and upper limits. Dashed lines show the effects of X-ray obscuration in the observed Lx/L^^m 
ratio. The thick solid, dot-dashed and dotted lines show the expected Lx / i-i2/jm for star-forming galaxies, 
considering different recipes to convert star formation rates into X-ray and IR luminosities. Right panel: 
Space density of CT AGN as a function of redshift. Measurements obtained by Tozzi et al. II15I1 . Alexander 
et al. lfl6ll . Fiore et al. [17] and this work are shown by squares, stars, pentagons and triangles respectively. 
Solid lines show the expected space density of CT AGN from the luminosity functions of Delia Ceca et al. 
ifisll . Red and blue symbols show measurements and expectations for Lx>l 0^^ and > 1 O^'* erg s ^ ' sources . 
These observations indicate a strong increase in the number of high-luminosity CT AGN at z>2. 
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